#ifndef _PO_MOM
#define _PO_MOM

#include <complex>
#include <blitz/array.h>
#include <blitz/tinyvec-et.h>
typedef std::complex<double> dComplex;
typedef blitz::TinyVector<double,3> Position3D;
typedef blitz::TinyVector<dComplex,3> Complex3D;
typedef blitz::Array<dComplex,2> Matrix2D;

extern blitz::Array<Position3D,1> vWire;
extern blitz::Array<double,1> Pattern;
extern blitz::Array<dComplex,1> V,I;
extern Matrix2D Z;
extern int iWire;
extern double DeltaLength;
extern double WireRadius;
extern const double PI;
extern const double Epsilon;
extern const double Mu;
extern const double C;
extern const double Freq;
extern const double Omega;
extern const double WaveNumber;
extern const double Eta;

#endif
